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ABSTRACT 

We compare the use of galaxy asymmetry and pair proximity for measuring galaxy merger fractions 
and rates for a volume limited sample of 3184 galaxies with —21 < Mb — blogh < —18 mag. and 
0.010 < z < 0.123 drawn from the Millennium Galaxy Catalogue. Our findings are that: 

(i) Galaxies in close pairs are generally more asymmetric than isolated galaxies and the degree of 
asymmetry increases for closer pairs. At least 35% of close pairs (with projected separation of less 
than 20 kpc and velocity difference of less than 500 km s~^) show significant asymmetry and are 
therefore likely to be physically bound. 

(ii) Among asymmetric galaxies, we find that at least 80% are cither interacting systems or merger 
remnants. However, a significant fraction of galaxies initially identified as asymmetric are contami- 
nated by nearby stars or are fragmented by the source extraction algorithm. Merger rates calculated 
via asymmetry indices need careful attention in order to remove the above sources of contamination, 
but are very reliable once this is carried out. 

(iii) Close pairs and asymmetries represent two complementary methods of measuring the merger rate. 
Galaxies in close pairs identify future mergers, occurring within the dynamical friction timescale, while 
asymmetries are sensitive to the immediate pre-merger phase and identify remnants. 

(iv) The merger fraction derived via the close pair fraction and asymmetries is about 2% for a merger 
rate of (5.2 ± 1.0) x lO"'' Mpc'^ Gyr"!. These results are marginally consistent with theoretical 
simulations (depending on the merger time-scale), but imply a flat evolution of the merger rate with 
redshift up to z ^ 1. 

Subject headings: galaxies: interactions - galaxies: structure 



1. INTRODUCTION 

In the A-dominated Cold Dark Matter (CDM) model, 
galaxies form via a process of hierarchical merging, in 
which more massive objects are assembled gradually via 
mergers of increasingly more massive subunits at pro- 
gressively lower redshifts (see Baugh 200fi for an intro- 
duction) . Observations of the merger rate of galaxies and 
its evolution with redshift provide important benchmarks 
for comparison with theoretical predictions, according to 
which massive galaxies have undergone several mergers 
between z ~ 3 and the present epoch, d oubling their 

' mass over the last half o f the Hubble time (jMurali et al.l 

' [200l lMaller^raII[200il . 

Unfortunately, it is generally difficult to ascertain ob- 
servationally if two galaxies are going to merge in the 
near future, or if and when a galaxy last experienced 
a merger. Galaxies undergoing mergers or interactions 
may exhibit peculiar morphological disturbances, such 
as tidal tails (|Arpll 19661 fToomre L ToomrdI 1972t iToomrd 
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|1977() . while merger remnants may also show relic struc- 
tures fr om their pre vious interactions, such as shells or 
ripples (|Malin fc C arter 1980, 1983). A more quantita- 
tive version of this approach measures the asymmetry 
of a galaxy's light distribution to identify interacting 
systems and merger remnants and estimate the galaxy 
merger rate, e.g. using t he CAS system (Bershad y, Jan- 
gren & Conselice 2000, IConselicg|2003t IConsclice et al.l 
|2003| and references therein) . 

Before galaxies merge they will be distinct systems 
with some radial velocity difference and projected sep- 
aration on the sky from which a merger time-scale can 
be calculated based on dynamical friction arguments. 
Galaxy pair statistics therefore could, in principle, pro- 
vide a powerful probe of mergers and their frequency. 
With modern imaging and spectroscopic surveys it is 
now possible to select pairs of galaxies which are not 
only close on the sky, but also at very similar redshifts, 
and are therefore more likely to be gravitationally bound. 
Such objects are referred to as 'dynamically close pairs'. 
The mathematical formalism to measure the close pair 
fraction and to deriv e the merger rate from it was devel- 
oped bv lPatton et a l. (2000) for the Second Southern Sky 
Redshift Survey and iPatton et al.l ()2002l ) for the Cana- 
dian Network for Observational Co smology 2 (CN0 C2) 
survey and later apphed to DEEP 2 (|Lin et al.ll2004D and 
the MGC (|De Propris et al.ll2005( l. Both of these latter 
studies found a low local merger rate, suggesting that 
only a small fraction of massive galaxies are formed by 
major mergers at z < 1. 

There are therefore two methods for observationally 
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measuring the merger history - finding galaxies in pairs 
that will eventually merge, and locating highly asymmet- 
ric galaxies. Neither of these two approaches is free from 
ambiguity however. For example, not all close galaxy 
pairs will be dynamically bound and result in mergers. 
The observational problem is that it is impossible to con- 
strain the projected relative motion of the two members 
of a close pair or to precisely dete r mine the radial dis- 
tance between them. iPatton et al.l (|200C0 estimate that 
close pair samples may suffer from contamination by as 
many as 50% unphysical pairs (superpositions). 

On the other hand, not all asymmetric galaxies are 
imminent or recent major mergers, as it is possible to 
induce significant asymmetries by other means (e.g., in- 
teractions, star formation, minor mergers, etc. - see, e.g., 
Conselice, Bershady & Jangren 2000a). Care must there- 
fore be taken when selecting mergers based on asymme- 
try (|Consehcdl2003l ). 

The two techniques - asymmetry and close pairs - 
are also likely to prob e different stages of a merger: 
iHernandez- Toledo et al.i (2005) suggest that galaxy mor- 
phology is very robust, and that significant asymmetries 
are induced only by a close interaction. Furthermore, 
N-body simulations indicate that a merger remnant may 
remain significantly asymmetric for an extended period 
after the merger event (Conselice 200_6,). Hence asymme- 
try is likely to probe the immediate pre-merger and an 
extended post-merger phase. In contrast, galaxy pairs 
will select all stages of the pre-merger phase except the 
very last, where the two galaxies can no longer be identi- 
fied unambiguously as individual objects in the imaging 
data. 

It is important to understand how galaxy pairs and 
asymmetric galaxies measure the merger process be- 
cause, for observational reasons, the two methods are 
preferentially employed in different redshift regimes. At 
high redshift, highly complete spectroscopic samples 
are difficult to obtain whereas relatively high resolution 
imaging is available from HST. Hence selection by asym- 
metry is often the m ethod of choice in this regime (e.g. 
IConselice et al.l [20031 ). At low redshift, the situation is 
essen tially reversed so th at the pair method is preferred 
(e.g. IPatton et al.ll2002t l. We clearly require a cross- 
calibration of the two methods on the same dataset be- 
fore drawing reliable conclusions regarding the evolution 
of the merger rate at high and low redshifts. 

In this paper we attempt to deliver this cross- 
calibration by searching for dynamically close pairs of 
galaxies and measuring the structural asymmetries on 
the same sample of objects. We then determine if paired 
galaxies are asymmetric, and also how the merger frac- 
tions and rates derived from both methods compare. A 
few questions we address include: are galaxies in close 
pairs more asymmetric than isolated galaxies? Does the 
degree of disturbance correlate with pair separation? Are 
more asymmetric galaxies more likely to lie in pairs, with 
the most asymmetric ones in closer systems? What is the 
distribution of galaxy asymmetries for the local universe? 
Are there isolated galaxies with high asymmetry, and are 
t hese merger remnan ts? 

IPatton et"al] (|2005| ) carried out a similar study by com- 
paring the asymmetries of isolated galaxies and those in 
close pairs, taken from the CN0C2 redshift survey, us- 
ing Hubble Space Telescope snapshots of close pair can- 



didates and a small sample of isolated CN0C2 galaxies 
falling within the HST field of view. For these galaxies, 
they measure asym metries using the Rt + Ra index of 
iSchade et all ([19951) . Galaxies in CN0C2 pairs are more 
asymmetric than isolated ones, but over only a limited 
range of pair separations (up to 21 k pc in projec- 
tion) se t by the original pair sel e ction by IPatton et al.l 
(|2002D . iHernandez-Toledo et al] (|2005 h have compared 
asymmetries in spiral-spiral pairs, ultraluminous infrared 
galaxies (co mmonly regarded as in teracting and/or merg- 
ing svstems: iLonsdale. F'arrah fc Srnith 20061 and a local 
sample of objects from lFrei et al.l () T996f ). showing that 
asymmetry increases with decreasing pair separation, al- 
though the samples were somewhat heterogeneous. 

To measure the merger rate from asymmetry and dy- 
namically close pairs requires a dataset with both high- 
quality imaging and very complete redshift i nformation. 
The Millennium Galaxy Catalogue (MGC; iLiske et al.l 
120031 ) is particularly well-suited to this task. The MGC 
consists of deep CCD imaging of a long, thin (36' x 
72°) equatorial strip, reaching to i?MGc = 24 mag and 
/is = 26 mag arcsec"^, and coinciding with both the 
2dF Galaxy Redsh ift Survey northern strip (2dFGRS; 
IColless et aLll200lD and with the Sloan Digital Sky Sur - 
vey data release 1 region fSDSS: 'Abazaii an et al.ll2003D . 
The MGC also includes a parallel redshift survey, reach- 
ing 99.8% complet eness at B^uic = 19 .2 mag and 96.0% 
at SjiGc: = 20 mag (jDriver et al.ll2005D . 

Because of its homogeneous high-quality imaging, 
its large, contiguous field of view, and high redshift 
completeness, the MGC enables one to compute the 
merger rate from the sta tistics of dynamically close pairs 
()De Propris et al.l l2005f ) and galaxy asymmetries. The 
main purpose of this paper is to carry out a comparison of 
these two approaches using the same, homogeneous, sam- 
ple of galaxies. As in previous MGC papers, we adopt 
a cosmology with Q,m = 0.3, IIa = 0.7 and calculate all 
distances with reference to Hq = 100 h km s~^ Mpc""'^. 

2. METHODOLOGY 
2.1. Sample 

In order to compare asymmetries and close pairs as 
indicators of the merger rate, we first select a volume- 
limited sample of galaxies from the MGC. This sam- 
ple includes 3237 galaxies with Bmgc < 20 mag, —21 < 
Mb - 51og/i < -18 mag and 0.01 < z < 0.123. These 
limits are shown in Figure[T] From this sample we extract 
close pairs with projected separations of up to 100 /i^^ 
kpc and isolated galaxies (i.e., with no companion within 
the specified magnitude and redshift limits and projected 
separation criteria). We also compute asymmetries for 
all galaxies in the volume-limited sample. Because of 
the absolute magnitude range considered, we can iden- 
tify mergers with a luminosity ratio of up to 1 : 16 for 
the brighter galaxies [Mb = —21 mag.), with the sam- 
ple becoming more incomplete for fainter galaxies. The 
range in luminosity ratios that can be explored biases the 
sample towards pairs with nearly equal luminosities. 

2.2. Galaxy Pairs 

We identify dynamically c lose pairs usi n g the above 
sam ple in the same wa y as IPatton et all pOOOl . I2002D 
and lPatton et al.l (|2005( ). However, we consider a wider 
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Fig. 1. — The volume-limited box used to select galaxies for 
analysis. The points are all MGC galaxies with Bugc < 20 mag. 



range of projected separations than iPatton et al.1 (|2005l ) 
because of our wider angular coverage and more complete 
spectroscopy. Two galaxies in our sample form a dynam- 
ically close pair if their projected separation is less than 
100 h^^ kpc, and their velocity difference is less than 
500 km s~^. We consider only the closest companion 
(multiplets are ignored) and galaxies with no such com- 
panion are regarded as isolated for our purposes. 

The above criteria for position and velocity separation 
were chosen as follows. The pairwi se velocity dispersio n 
of 2dF galaxies is ~ 500 km s'^ (jHawkins et alll2003f ): 
therefore objects with larger velocity separ ations are un- 
likely to be bound (see also Figure 4 in iPatton et al.l 
|2000( ). In addition, simulations show that pairs with 
larger velocity separations are unlikely to merge (Carl- 
berg, Pritchet & Infante 1994), while a velocity sepa- 
ration of 500 km s~^ is generally larger than the in- 
ternal velocity dispersion of the galaxies under study, 
which decreases the likelihood of a significant interaction 
(jMakino & Hut 1997). 

Simulations also show that galaxies with similar 
masses, and with a projected separation of le ss than 
20 kpc are likely to merge within 1 Gyr (jBarnesI 
[19881 iCarlberg et al.lll99l lConsehcel[200l . The upper 
limit on the projected separation we adopt in this paper 
(100 kpc) is arbitrary, but it should isolate all ob- 
jects likely to merge in the near future, as well as those 
wh ere the interaction has induced significant asymmetry 
(of. iHernandez-Toledo et al.ll2005| ). 

It is of course possible that we miss some companions 
due to our incomplete redshift coverage. There are 399 
Bmgc < 20 mag galaxies (4% of all B^cc < 20 systems) 
in the MGC for which we have not obtained a redshift. 
We use the projected separation criterion to calculate 
whether any of these galaxies could be the closest com- 
panion to a galaxy in the volume-limited sample. If so, 
we assume that the galaxy without redshift has the same 
redshift and fc-correction as its putative companion, and 
require that it obeys the absolute magnitude cuts to be 
part of the volume-limited sample. Objects in the sam- 
ple for which the closest companion may be one of the 
399 galaxies without redshifts are however not included 
in our analysis. This decreases our sample to 3184 galax- 
ies from the original 3237. This contains 112 galaxies in 
pairs with projected separation of less than 20 h^^ kpc, 
and 2561 isolated galaxies. While all pairs within the 
specified luminosity and redshift limits are included, the 
sample is incomplete towards pairs with large luminosity 



differences and tends to be more complete for pairs of 
nearly equal luminosities. The median magnitude differ- 
ence for the above 112 pairs is 0.6 mag., equivalent to a 
typical 1:2 luminosity ratio (i.e., we primarily measure 
major mergers). 

2.3. Galaxy Asymmetries 

Asymmetries for galaxies within the volume-limited 
sa mple are calc u lated using both the Rt + Ra index 
of ISchade et al.l (I1995D and th e more recent CAS in- 
dices of IConselice et alj ()2000af) . which have been cali- 
brated on local galaxy samples to identify merger rem- 
nants. The details on how to calcula te the asym- 
metrie s us ed in this paper are given in ISimard et al.l 
(I2002D and I Allen et all (l2006f) . f or the Rt + i^A method , 
and IC onselice et al.' ('2000a'); 'Bershadv^eTaD (j2000D : 
IConsel ice (2003); Consclicc et al. (2004) for the revised 
CAS indices. The CAS (concentration, asymmetry, 
dumpiness) parameters are a non-parametric method for 
measuring the structures of galaxies resolved on CCD 
images. One particular feature of this system is that 
galaxies with extremely high asymmetries, usually A > 
0.35, are likely to be in a major (post-)merger phase. 
This has been calibrated with asymmetry measurements 
of nearby and distant normal and merging galaxies 
jConsehce et al."2000a"; |Consehcel [2003t IConsehce eTel] 
|2005; Bridge et al. 200 3), ancTof 'galax ies' seen in N- 
body merger simulations ('Consel ic3l2006f l . Furthermore, 
Conselice, Bershady & Gallagher (2000b) have shown 
that galaxy mergers have both high asymmetries and 
high HI line width asymmetries. 

Our method for computing the CAS parameters is 
slightly different from that in previous works. Because 
we are examining nearby bright galaxies in a general 
field survey, our sample was not specifically selected 
to be clean of foreground and background contamina- 
tion. Our sample therefore represents a more gener- 
alised case for measuring asymmetries in field surveys 
of nearby galaxies. For e ach galaxy we use a SExtractor 
(|Bertin fc Arnouti I1996D segmentation map to replace 
all nearby objects with the sky, using the same noise 
properties as measured directly from the background. In 
this way we eliminate the possible contamination of the 
asymmetry measurement by these objects. This method 
is however not always perfect when there are very nearby 
bright stars or large galaxies. We must therefore keep in 
mind that any sample of highly asymmetric galaxies may 
be contaminated by galaxies that have an intrinsically 
low asymmetry, but whose asymmetry measurement was 
corrupted by a nearby object. 



3. RESULTS 



3.1. 



The Asymmetry Distribution 

F igure|2|shows the distribution of asymmetries for both 
the lSchade et al.l ([l995) and CAS methods: we plot the 
fraction on the ordinate logarithmically, in order to em- 
phasise the small number of more asymmetric galaxies. 

For the Rt + Ra index we find that the asymme- 
try distribution peaks at ~ 0.05, with a long tail of 
more asymmetric o bjects. Following the definitions of 
IPatton et all (|2005l ). (18.8 ± 0.8)% of galaxies are 'asym- 
metric' (i?T + Ra> 0.13) and (12.4 ± 0.7)% are 'highly 
asymmetric' {Rt+Ra > 0.16). However, unlike the CAS 
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Fig. 2. — The distribution of asymmetry values for all galaxies 
in the volume-limited sample. Note that the ordinate is plotted 
logarithmically to better show the small number of asymmetric 
galaxies. 

indices, Rt + Ra has not been calibrated to specifically 
identify merging systems or merger remnants. Visual 
examination of the 2% most asymmetric objects (by the 
Rt + Ra criterion) shows that about 1/3 are in close 
pairs, 1/3 seem to be merger remnants and the rest ap- 
pear to be edge-on galaxies or classical irregulars and 
late-type spirals (see also section below). 

Among the close pairs with projected separation of less 
than 20 kpc we find that (25.5±5.5)% are 'asymmet- 
ric' and (16.9 ±4.3)% are 'highly asymmetric'. The pair 
fractions and merging pair fractions (assuming that 50% 
of the pairs actually merge ) are in good agreement with 
the values derived bv .Patton et al.l (2005I ). but we find a 
much larger fraction of asymmetric or highly asymmet- 
ric isolated galaxies. This is likely due to the re l ativel y 
small sample of isolated galaxies in iPatton et al.l (|2005l ) , 
the lower physical resolution of the comparatively shal- 
low HST imaging used (leading to a smoothing of the 
low surface brightness brightness asymmetric features), 
different isolation criteria used and treatment of segmen- 
tation maps. 

The CAS method defines galaxies with A > 0.35 as 
unusually asymmetric: Conselice et al. (2000a) suggest 
that such objects may be interpreted as merger rem- 
nants. The distribution of CAS A values in Figure [2] 
peaks at ~ 0.11, suggesting that A is sensitive to rel- 
atively minor deviations from symmetry. We find that 
(4.1 ± 0.4)% of our sample have A > 0.35 with this bhnd 
selection. At face value, and assuming that 50% of the 
close pairs (with separation less than 20 kpc) merge, 
this is about a factor of two larger than the 'merging 
pair' fraction. However, before we draw any conclusions 
regarding possible discrepancies between the two meth- 
ods we will need to evaluate more carefully what types 
of objects have been selected by the A > 0.35 cut. 

Figure [3] compares asymmetry measurements for galax- 
ies in various isolation classes (from close pairs to galaxies 
we regard as being isolated) for both CAS A and Rt+Ra 
and shows the limits used to distinguish asymmetric or 
highly asymmetric galaxies for both methods. Not sur- 
prisingly, Rt + Ra thresholds appears to be considerably 
less stringent than CAS A, although there is a broad re- 
lation between the two indices. One can see that there 
are several galaxies in close pairs with separations less 
than 40 kpc above the A = 0.35 line, as well as 
numerous isolated galaxies, while only a few pairs with 
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Fig. 3. — Comparison of CAS A values and Rt + Ra for galaxies 
in different separation classes. The horizontal thick dashed line 
indicates the A = 0.35 asymmetry limit in the CAS system. The 
two vertical thick dashed lines show the Rt + Ra = 0.13 and 0.16 
limits. 

larger separation have high CAS asymmetry. There are 
numerous objects deemed as highly asymmetric by the 
Rt + Ra index, but with A < 0.35. In most cases, these 
are edge-on or late-type spirals which have not been de- 
contaminated from the sample as we do for the high A 
systems. 

3.2. The Nature of Asymmetric Galaxies 

A blind study such as this requires that we examine by 
eye those systems which are asymmetric. There are two 
reasons why one should not simply assume that all MGC 
galaxies with A > 0.35 are imminent mergers or merger 
remnants. First of all, as pointed out in Section [^751 the 
sample may contain objects whose asymmetry measure- 
ments have been corrupted by nea rby br ight /large ob- 
jects. Secondly, the A value of 0.35 (|Consel ice 2003) was 
established using the sample of normal nearby galaxies of 
iFrei et all ()1996( ) and a collection of starburst galaxies, 
which were selected to span a large range of morphologi- 
cal types. In contrast, here we are using a volume- limited 
sample. Also, although the CAS A is reasonably robust, 
it does depend on the physical resolut ion of the data 
(|Conselice et al.l l2000at IConselicd l2003l ) . Therefore we 
have visually examined all 129 galaxies with A > 0.35 in 
our sample and classified them into the following cate- 
gories: 

• Likely merger remnants: although this is a 
somewhat subjective classification, these objects 
appear to be relatively obvious 'trainwrecks' or 
very disturbed disks with multiple, similarly bright 
light components. We find 23 objects in this cate- 
gory, and we show MGC postage stamp images of 
these galaxies in Figure [H 

• Likely imminent mergers: these are disturbed 
galaxies similar in appearance to the above, but 
satisfying two additional requirements: (i) they are 
a member of a close pair as identified above, i.e. 
have a 'major' companion; (ii) they show evidence 
of interaction with their companion. We find 13 
such objects and show their images in Figure [51 
The projected separation to the companion, d, in 
this sample ranges from 2.5 to 34 h^^ kpc. 
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Fig. 4. — Postage stamp MGC images of galaxies with A > 0.35 which are regarded as merger remnants. The scale of each image is 
indicated by the horizontal bar which is 10 arcsec. 



• Other: there are 8 galaxies that are not immi- 
nent mergers and whose asymmetry measurements 
are not contaminated (see below). However, 4 of 
these galaxies are faint and/or low surface bright- 
ness objects which we cannot reliably classify, and 
all of these look unusual and could be merger rem- 
nants. In contrast, the other 4 galaxies in this cat- 
egory are most likely not recent merger remnants: 
one is a face-on disk, one is a spheroidal galaxy 
(with some low-level debris nearby), and the other 
two are minor mergers (as determined from the full 
MGC catalogue). 

• Contaminated: galaxies where the asymmetry 
appears to be due to the presence of a bright star 
in the vicinity. 28 galaxies are thus affected. Fur- 
thermore, for 57 galaxies SExtractctor did not con- 
struct correct segmentation maps, which led to er- 
roneously high CAS A values. If these objects are 
(wrongly) included in the sample of highly asym- 
metric galaxies they account for 50% of the sam- 



ple. Hence, without paying particular attention to 
the construction of the segmentation maps, and/or 
visual vetting of the sample, a large fraction of 
A > 0.35 galaxies turn out to be false positive 
mergers. 

If we exclude objects with contaminated measurements 
and with bad segmentation maps, we find that 80% 
of the asymmetric galaxies in our sample are obvious 
imminent or recent mergers. The status of the remainder 
is more uncertain. However of the remaining 8 systems, 
four are possible merger remnants. 

We note that over half of all systems initially found 
with A > 0.35 were contaminated by stars, or were over- 
segmented. This demonstrates that vetting of asym- 
metry selected merger samples is advisable (see also 
IConselic^l2003[ ) . Although contamination by stars is less 
of a problem at high rcdshifts, where, presumably, the 
asymmetry index will be applied in future surveys, con- 
tamination b y neighbouring unrelat ed galaxies may be 
an issue (e.g. iKampczvk et al.|[2007D . 
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Fig. 6. — As Figure|4]but for galaxies with A > 0.35 which cannot be reliably classified. 



As noted above, we found from our visual inspec- 
tions that 57 of our initial sample of A > 0.35 galax- 
ies were asymmetric due to Sextractor having produced 
bad segmentation maps for these objects. Most of these 
were late- type, face-on or edge-on spirals with prominent 
star- forming regions and/or spiral arms, and SExtractor 
falsely separated these into multiple sources. As noted 
in Section 2.3, the CAS code attempts to avoid the cor- 
ruption of a galaxy's CAS parameters through nearby 
objects by replacing these with blank sky. In the case 
of the shredded galaxies this procedure created 'holes' in 
their light distribution and hence artificially high asym- 
metries. To deal with this we re-ran the CAS code on 
all galaxies without using any segmentation maps at all. 
Each galaxy was then assigned the smaller A value from 
the two runs, and 57 galaxies went from A > 0.35 to 
A < 0.35. 



Hence, in order to avoid significant contamination of 
high-asymmetry samples with late-type disks in future 
large surveys (where visual inspection will be unfeasible) , 
care must be taken when constructing the segmentation 
maps and/or when masking neighbouring objects. 

3.3. Asymmetry and Pair Separation 

Figure [7] shows the values of Rt + Ra and CAS A 
for galaxies in pairs, plotted against the projected sep- 
aration of the pair members. Galaxies in pairs with 
c? < 50 kpc clearly show an excess asymmetry over 
galaxies in pairs with wi der separations. This appears to 
conf irm the findings of iHernandez-Toledo et al.l (|2005l ) 
and iPatton et al.l ()2005f ) that being in a close pair in- 
duces asymmetries (via tidal stresses and star forma- 
tion). Hence there is no doubt that at least a fraction 
of the close pairs consist of physically bound objects un- 
dergoing an interaction. 
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Fig. 7. — Galaxy asymmetry versus pair separation using both 
the Rt + Ra a-nd the CAS A indices. 
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Fig. 8. — Fraction of asymmetric galaxies in five projected sep- 
aration bins. The horizontal lines show the fraction of isolated 
galaxies with i?y + Ra > 0.13 (dashed line) or > 0.16 (dot-dashed 
line) or A > 0.35 (dashed line in upper panel). 

In order to show this more clearly we plot the frac- 
tion of 'asymmetric' galaxies vs. pair separation for five 
20 kpc wide bins of projected separation in Figure[8l 
We also indicate the fraction of asymmetric galaxies (ac- 
cording to both methods) among the isolated galaxies. 
The fraction of asymmetric galaxies increases with de- 
creasing pair separation, consistent with the interpreta- 
tion that closer pairs result in a progressively stronger 
interaction. 

We also compare the cumulative distribution of asym- 
metries for the five projected separation bins used above, 
and for isolated galaxies in Figure [Q] As implied by Fig- 
ures [7] and [51 we see an excess of asymmetric galaxies over 
the distribution for isolated galaxies only in the two low- 
est separation bins. However, in the smallest separation 
bin, there is an excess fraction of asymmetric galaxies 
even at the highest asymmetries. This implies that, at 
least in the closer pairs, we are witnessing a significant 
degree of morphological disturbance. 

While these figures show that Rt + Ra is also capable 
of identifying mergers and merger remnants, this index 
has not been calibrated with the same thoroughness as 
the CAS index, where extensive work has been carried 
out by Conselice and collaborators. For this reason, from 
now on, we use the CAS index exclusively to derive the 
merger rate and compare to the pair fractions. 



By comparing the differential CAS A distribution of 
isolated galaxies with that of galaxies in close pairs with 
d < 20 kpc we can derive a limit on the fraction, 
/, of unphysical pairs in this sample. For our purposes, 
we assume (or define) galaxies in 'unphysical' pairs to 
be those that have the same A distribution as isolated 
galaxies. This calculation will give us therefore a lower 
limit on the fraction of paired galaxies that are in in- 
teracting systems. We denote the normalised differential 
A distribution of isolated galaxies by i{A), and the un- 
known distribution of galaxies in true, physical pairs as 
t{A). Hence the normalised A distribution of all appar- 
ent galaxy pairs with d < 20 kpc can be written 
as: p{A) = f i{A) + (1 - /) t{A). Since / must he in 
the range < / < 1 and since t{A) > for all A, we 
find that p{A) > f i{A) must hold for all A. By com- 
paring the observed distributions p{A) and i{A) we find 
that the maximum value of / for which this condition 
still holds is 0.65, where this constraint is driven by the 
low- A end of the distributions {A < 0.15) as expected. 
In other words, of all galaxies in dynamically close pairs 
with d < 20 kpc at least 35% lie in physical pairs. 
Th is compares well wi th the fraction of 50% estimated 
by iPatton et al.l (|2000f) . and we will use this number in 
the following calculations. 

4. PAIR AND ASYMMETRY MERGER TIME-SGALES: 
IMPLIGATIONS FOR THE MERGER FRAGTION AND 
GALAXY EVOLUTION 

We are now in a position to calculate the merger time- 
scale and fraction for close pairs and asymmetry-selected 
galaxies, and compare the lifetime of the merger remnant 
to the merger time-scale for the progenitor pair. 

In Section FS.SI we found that only pairs with small sepa- 
ration {d < 50 h^^ kpc) are likely to merge within a short 
time interval, and that among these the majority of likely 
mergers are found among pairs with d < 20 h^^ kpc 
(based on the presence of asymmetries as an indica- 
tor of an on-going interaction). To a first approxima- 
tion the close pair fraction for these latter objects can 
give a measure of the merger rate. We find that the 
pair fraction is (3.5 ± 0.4)%. In order to correct for 
pairs which are missed because of redshift incomplete- 
ness we ad opt the 18.5% incomple teness correction cal- 
culated by iDe Propris et al.l (|2005| ) for MGC close pairs 
with d < 20 kpc, yielding a close pair fraction of 
(4.1 ± 0.4)%. Assuming that 50% of these close pairs are 
likely to be physical pairs and hence future mergers, we 
derive a merger fraction of (2.1 ± 0.2)% for our volume- 
limited sample. 

In Section [3.21 we found that high asymmetry identifies 
both merger remnants and to a lesser extent future merg- 
ers. To derive a merger fraction that is directly compara- 
ble to the one derived from the close pairs above (which is 
essentially a 'progenitor galaxy' fraction), we must take 
into account that each merger remnant was produced 
by (at least) two progenitor galaxies and hence count it 
twice. In addition we apply the above incompleteness 
correction to the 8 objects in pairs with d < 20 kpc. 
This yields a galaxy merger fraction of(1.9±0.2)%. 

The similarity between the merger fractions derived by 
these two methods suggests that the total time-scale for 
asymmetry (i.e. including the pre-merger phase) is sim- 
ilar to the merger time-scale of a d < 20 h"^ kpc pair. 
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shows the observed asymmetry distribution for isolated galaxies (shown in the bottom right-hand panel). 



We can determine the pre- and post-merger fraction of 
the total asymmetry time by comparing the number of 
imminent mergers (corrected for redshift incompleteness) 
with the number of merger remnants (x2) identified by 
asymmetry. We deduce that the pre-merger phase takes 
up 24% of the total asymmetry time, consistent with re- 
sults derived from N-body simulations (Conselice 2006). 
We notice, however, that there are several galaxies in 
pairs with A > 0.35 that have large separations (some as 
large as 34 h^^ kpc). It is likely that these objects are 
on parabolic or hyperbolic orbits or may not merge until 
after several more passages. 



The general picture that emerges from these results is 
the following: close pair selection (with d < 20 kpc) 
is sensitive to most stages of the 'pre-merger'. From 
dynamical friction considerations (see equation 7 in 
[ConscUcc 2006) we find that the merger time-scale for 
a 'typical' close pair in our sample, with average sepa- 
ration of 12 kpc and average velocity separation of 
150 km s^^ is about 0.3 Gyr. Being in a close (physical) 
pair enhances asymmetry and during the last ~ 22% (or 
65 Myr) of this pre-merger phase the asymmetry has be- 
come so significant that the imminent merger is selected 
by the A > 0.35 cut. The total asymmetry time is much 
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longer than this because the remnant remains asymmet- 
ric for 0.21 Gyr after the two systems inerge, and the new 
system dynamically relaxes (jConselicd I2006D . Hence, as 
implied by N-body simulations, asymmetry identifies the 
immediate pre-merger stage as well as an extended post- 
merger phase. While the absolute time-scales above are 
quite uncertain their ratios are not, but are fixed by our 
results, i.e., if the merger time-scale is revised to 1 Gyr 
then the asymmetry time-scales must also change accord- 
ingly. 

A comparison with theoretical predictions is somewhat 
difficult, as there may not be a one-to-one correspon- 
denc e between galaxies and their host dar k haloes (see, 
e.g.jBerrier et al.ll2006l ). With this caveat. iMaller et all 
(|2006D predict a merger rate of 0.054 mergers per Gyr 
for massive galaxies with mass ratio 1:1 to 1:3. From 
our data we find that 66% of the MGC pairs consist of 
galaxies within this mass ratio, yielding a pair fraction 
of ~ 0.014. The merger rate we derive is marginally con- 
sistent with the above prediction for a merger time-scale 
of 0.3 Gyr or too low by a factor of 3.9 for a time-scale 
of 1 Gyr, similar to the comparison shown by Bell et al. 
(2006 - their figure 2) from theoretical models. 

The merger r ate we derive, following equation 3 of 
iLin et all (f2004), is (5.2 ± 1.0) x 10"* Mpc'^ Gyr"!. 
Strictly speaking this is a lower limit, as some pairs 
are certainly missed, but it should be a fair estimate 
of the major merger rate for nearly equal luminos- 
ity pairs, which are less incomplete. For compari- 
son, [LiZ^iE^n l|2004[ ) derive an average merger rate of 
4 X lO"'' Mpc"^ Gyr"^ at 0.5 < z < 1.2 for galax- 
ies with —21 < Mb < —19, which is abo ut one or- 
der of magnitude lower than that derived by IConselic^ 
(|2003f) . but in good agreement with our z ~ 0.1 value. 
It is much larger, by almost two order of magnitude, 
than the merger rate be tween luminou s red g alaxies at 
z < 0.36 measured by iMasiedi et al.l ()2006[ ). suggest- 
ing that the majority of mergers betw een galaxies i n 
the z < 0.5 universe are not 'dry' (c.f. iBlantonI l200l . 
Our results are consistent with a flat evolution of the 
merger rate, i.e., as (1 + z)~^ ()Lin et al.l [200l rather 
than the (l-|-z)^ dependence expected from CDM models 
(iGottlober. Klvpin fc Kravtsovll200ll but see discussion 
bv iBerrier et al.ll2006[ ). If the merger rate evolution is 
flat out to z ~ 1 only a small fraction of massive galaxies 
can have formed via major mergers in the last 1/2 of the 
Hubble time. 

5. SUMMARY 

In this paper we have compared the close pair fraction 
and the fraction of asymmetric galaxies in a volume- 
limited sample of galaxies drawn from the MGC. Our 
sample consists of 3184 galaxies between —21 < Mb — 
\ogh < —18 mag and 0.01 < z < 0.123. The main con- 
clusions of this work are: 

• The CAS selection for mergers, A > 0.35, after cor- 
rection for contamination, is highly successful at lo- 
cating galaxies in major mergers with a corrected 
fraction of ^ 80%. The use of a blind selection how- 
ever can result in many false positives. To create a 



clean merger sample selected with the asymmetry 
index it is necessary to either visually examine all 
high-A galaxies, or use othe r automated me thods, 
such as a high dumpiness (jConselicd l2003l) , and 
axis ratios to remove contamination. 

• We find a connection between pairs and asymme- 
try. Galaxies in pairs tend to be more asymmetric 
and the asymmetry tends to increase with decreas- 
ing pair separation, implying that these galaxies 
are actually interacting. We estimate that at the 
very least 35% of galaxies in close pairs with pro- 
jected separation of less than 20 kpc are actu- 
ally physically bound based on their excess distor- 
tion compared to isolated galaxies. 

• The overlap between our close pair and asymmetry- 
selected merger samples is small. Of all galax- 
ies in pairs with projected separation less than 
20 kpc only 7% were identified as likely merger 
events by the asymmetry method. Hence it is likely 
that the pair and asymmetry methods are tracing 
somewhat different phases and perhaps mass ratios 
of merging galaxies, while some objects in pairs are 
not asymmetric enough to be tagged by the CAS 
indices. 

• We find that after correcting the close pair sam- 
ple for chance superpositions and the asymmetry 
sample for contamination (by stars, and most im- 
portantly bad segmentation maps) that the merger 
fraction and rate for the two methods are very sim- 
ilar. This suggest that both methods are tracing 
the same underlying merger process, although at 
somewhat different phases of the merger process. 

• Once we have accounted for contamination and in- 
completeness in both methods we derive a consis- 
tent merger fraction of 1.9% per 0.3 Gyr and a 
merger rate of (5.2 ± 1.0) x lO"'' Mpc'^ Gyr"!. 
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